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Statement  of  Work 


1.  Determine  experimentally  the  conditions  under  which  plasmas  can  be  initiated  and 
sustained  in  the  three  energy  addition  modes  (TMoj,  TEq],  and  planar).  Initial 
testing  will  be  with  nitrogen  and  helium.  Nitrogen  will  be  used  to  simulate  the 
molecular  nature  of  hydrogen  while  helium  will  simulate  atomic  hydrogen.  Final 
testing  will  be  with  hydrogen.  Parameters  to  be  examined  include  gas  composition, 
pressures  and  flow  rates  and  microwave  power.  Quantify  system  heat  losses. 

2.  Measure  spectroscopically  electron  and  ion  temperatures  and  densities  in  microwave 
generated  plasmas  in  the  three  energy  addition  modes.  Due  to  high  gas  pressures, 
local  thermodynamics  equilibrium  (LTE)  will  be  initially  assumed  although 
nonequilibrium  effects  will  be  examined.  The  principal  measurement  will  be  ion 
temperature  as  this  translates  to  thrust.  The  effects  of  various  hot  and  cold  gas 
mixing  schemes  on  final  temperature  will  be  studied. 
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Summary  of  Results 


Free-floating  spherical  plasmas  have  been  generated  in  stationary  and  flowing 
nitrogen  and  helium  gas  contained  inside  a  10.2  cm  diameter  spherical  quartz  vessel 
located  within  a  cylindrical  resonant  cavity  operated  in  the  TMoj2  niode  at  a  frequency 
of  approximately  2.45  GHz^.  The  plasmas  are  approximately  two  inches  in  diameter  and 
are  centered  within  the  quartz  sphere.  Power  coupled  to  the  plasma  was  measured  as  a 
function  of  gas  composition,  flow  rate  and  pressure  and  microwave  power  input  to  the 
cavity  only  for  those  conditions  which  resulted  in  the  plasma  being  stabilized  in  the 
center  of  the  quartz  sphere  away  from  adjacent  walls. 

The  coupling  efficiencies  of  helium  and  nitrogen  plasmas  were  measured  with  no  gas 
flow  and  an  input  power  of  400  W.  Similar  measurements  were  made  for  a  helium  plasma 
with  a  gas  flow  rate  of  1.056  x  10"^  kg/s  and  400  W  input  power.  Figure  1  shows  the 
coupling  efficiency  versus  pressure  curves  for  the  three  cases  considered.  The  coupling 
efficiency  is  defined  to  be  the  percentage  of  input  microwave  power  actually  absorbed  in 
the  gas.  The  most  striking  feature  is  the  much  greater  range  of  operation  of  the  helium 
discharges.  The  nitrogen  discharge  is  limited  to  low  pressure  operation  (less  than  30 
kPA  (absolute)),  whereas  the  helium  discharges  operate  at  well  above  250  kPa  (absolute). 
The  helium  discharges  generally  exhibit  higher  coupling  efficiencies,  with  peak  values  of 
51%  (no  flow)  and  63%  (with  gas  flow).  In  comparison,  the  nitrogen  discharge  exhibits  a 
peak  coupling  efficiency  of  40%. 

For  a  given  plasma  in  the  resonant  cavity,  there  is  a  maximum  operating  pressure, 
above  which  the  plasma  is  extinguished.  Tests  were  performed  to  investigate  this 
stability  boundary  as  a  function  of  input  power.  The  first  tests  examined  nitrogen  and 
helium  plasmas  with  no  gas  flow,  with  later  tests  examining  a  helium  plasma  with  various 
gas  flow  rates.  As  is  shown  in  Figure  2,  the  maximum  operating  pressure  for  a  given 
input  power  is  significantly  greater  for  helium  discharges  than  for  nitrogen  discharges. 
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It  was  possible  to  sustain  a  nitrogen  discharge  up  to  a  pressure  of  only  40  kPa  (absolute) 
with  an  input  power  of  450  W,  whereas  helium  discharges  displayed  much  higher  maximum 
pressures  for  cases  both  with  and  without  flow  present. 

The  helium  discharges  for  various  flow  rates  all  have  similar  stability  boundaries, 
which  are  approximately  linear  at  first,  then  reaching  a  maximum  possible  operating 
pressure.  For  the  no-flow  case,  the  maximum  pressure  is  284  kPa  (absolute)  at  an  input 
power  of  400  W.  This  input  power  is  a  maximum  for  the  discharge  with  no  flow,  as  it 
was  found  that  further  increase  in  the  input  power  resulted  in  an  excursion  of  the 
discharge  to  the  surface  of  the  quartz  sphere. 

The  introduction  of  a  gas  flow  did  not  affect  the  stability  boundary  dramatically, 
but  did  raise  the  maximum  operating  pressure  of  the  discharge.  For  a  mass  flow  rate  1.5 
X  10”^  kg/s,  the  maximum  operating  pressur  was  318  kPa  (absolute),  at  an  input  power  of 
433  W.  The  tests  were  discontinued  at  an  input  power  of  505  W,  as  higher  input  powers 
resulted  in  the  plasma  positioned  close  to  the  quartz  surface.  For  an  increased  mass 
flow  of  1.056  x  lO"'^  kg/s,  the  maximum  operating  pressure  was  330  kPa  (absolute)  at  an 
input  power  of  400  W,  with  a  similar  profile  as  before. 

A  numerical  model  of  the  one-dimensional  planar  propagating  microwave  plasma  in 
hydrogen,  helium  and  nitrogen  gas  has  been  successfully  formulated. This  model 
numerically  integrates  the  system  of  governing  equations  consisting  of  the  one¬ 
dimensional  steady  energy  equation  and  Maxwell’s  equation  describing  the  propagation  of 
the  microwave  energy.  Due  to  thermal  conduction  of  the  cold  gas  ahead  of  it,  the 
plasma  propagates  toward  the  microwave  energy  source  at  a  velocity  determined  by  the 
energy  balance  between  the  absorbed  microwave  power  and  the  heated  gas  which  is 
convected  away  downstream. 

The  two  governing  equations  were  numerically  integrated  using  a  fifth/sixth  order 
variable  step  Runge-Kutta  scheme.  An  iterative  method  was  used  to  determine  the 


propagation  velocity  eigenvalue,  pu,  similar  to  the  method  used  by  Kemp  and  Root^  to 
solve  for  the  propagation  velocity  of  a  laser  heated  plasma.  The  propagation  velocity, 
maximum  temperature  and  percent  power  absorbed  were  calculated  as  functions  of  the 
input  microwave  power.  The  propagation  velocity  was  found  to  rise  with  increased 
microwave  power  but  the  maximum  gas  temperature  is  constant  for  hydrogen  and  nitrogen 
because  significant  dissociation  is  occurring  in  this  temperature  range  and  is  absorbing 
the  additional  absorbed  power.  The  maximum  gas  temperature  for  helium,  which  is  not 
dissociating,  was  found  to  rise  with  increasing  input  power.  It  was  found  that  all  of  the 
input  power  was  either  reflected  or  absorbed  with  no  power  being  transmitted  through 
the  plasma  and  that  the  present  power  absorbed  for  hydrogen  decreases  with  increased 
input  power  while  the  percent  power  absorbed  for  helium  and  nitrogen  remained  fairly 
constant  as  a  function  of  input  microwave  power.  Radiative  heat  loss  from  the  plasma 
was  subsequently  added  to  the  model  but  was  found  to  have  no  effect  on  the  numerical 
results.'^ 

The  results  of  the  experiment  are  shown  in  Fig.  3^  which  gives  the  measured  plasma 
velocities  for  helium  and  nitrogen  at  1  atm  together  with  the  numerical  predictions.  For 
helium  a  reasonable  agreement  between  numerical  and  experimental  data  could  only  be 
found  towards  lower  power  levels,  however,  an  enormous  deviation  can  be  observed  for 
power  levels  above  1550  W.  It  is  believed  that  in  this  power  range  a  change  in  the 
propagation  mode  takes  place  towards  resonant  radiation.  Nitrogen  doesn’t  show  this 
kind  of  behavior  and  agrees  reasonably  well  with  the  numerical  model.  The  percentage 
absorbed  power  values  for  helium  at  1  atm  range  from  the  mid  60  to  the  low  70  percent 
and  for  nitrogen  from  the  low  to  mid  60  percent.  These  values  were  much  higher  than 
those  predicted  by  the  model. 
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I.  STATtJS  OF  RESEARCH  EFTORT 

Our  research  errphasis  during  the  past  year  has  been  fcxaassed  on  two 
radiation-gasdyna]id.c  prcialems,  namely  solar  thermal  prcpulsion  and 
micrcwave  thermal  prcpulsion.  Primary  enphasis  was  on  the  solar  concept 
with  attention  being  shifted  tcward  the  microwave  concept  near  the  end  of 
the  contracted  effort.  The  purposes  in  both  areas  was  to  assess  the 
scientific  feasibility  of  the  particular  propulsion  conc^jt  on  the  basis  of 
detailed  analytical  models.  With  respect  to  the  solar  probloii,  the  key 
issue  is  to  identify  vAiether  or  not  fliaids  can  be  obtained  that  will  absorb 
the  solar  radiation  directly.  An  assessment  of  this  feasibility  is 
summarized  herein  and  presented  in  more  detail  in  published  articles.  The 
primary  issue  with  microwave  prt^xilsion  is  to  determine  vhether  or  not  the 
location  and  size  of  the  heated  plasma  region  can  be  controlled  adequately 
to  prevent  arcing  to  the  wall.  The  analytical  portion  of  this  study  is 
just  getting  underway  and  is  being  performed  in  conjunction  with  the 
ejq^eriroental  program  described  in  the  earlier  part  of  this  r^»rt. 

A.  Feasibility  of  Solar  Proculsion 

Solar  thermal  heating  appears  to  be  a  viable  and  realizable  prcpulsion 

conceit  if  the  solar  energy  can  be  coupled  directly  to  the  working  flxoid. 

This  is  referred  to  as  "direct”  solar  propulsion.  An  alternative  is  to  use 

the  solar  energy  to  heat  a  working  surface  of  the  engine  and  then  to  heat 

(1-3) 

the  working  fluid  indirectly  by  using  it  to  cool  this  surfaces  .  This 
"indirect"  solar  thermal  conc^^t  is  clearly  inferior  to  the  direct  concept. 
In  the  indirect  conc^t,  the  maximum  tenperature  within  the  cycle  occurs  on 
a  material  surface.  The  peak  tenperature  of  the  working  fluid  is  limited 
to  tenperatures  that  are  less  than  this  material  tenperature.  In  the 
direct  solar  thermal  concept,  the  maximum  tenperature  occurs  in  the  fluid 
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and  standard  techniques  can  be  used  for  cooling  adjacent  surfaces  to  keep 

them  at  Icwer  tenperatures.  Ihis  inplies  that  hi^er  working  temperatures 

can  be  used  resulting  in  inproved  thermodynamic  efficiencies.  In 

particular,  hi^er  tenperatures  in  propulsion  systems  correspond  to  hi^er 

specific  inpulses.  Ihe  inaximum  tenperature  in  direct  solcur  themal 

propulsion  thus  beccmes  limited  by  the  thermodynamic  characteristics  of  the 

solar  radiation,  \»hereas  in  the  indirect  case  it  is  limited  by  materieil 

considerations.  As  indicated  above,  the  major  question  concerning  the 

feasibility  of  solar  thermal  propulsion  is  vhether  or  not  the  solar  energy 

can  be  d^xDsited  directly  in  the  working  fluid.  It  is  toward  this  issue 

that  the  present  research  is  directed. 

Ihe  first  issue  concerning  the  feasibility  of  absorbing  solar  energy 

in  a  gaseous  working  fluid  is  to  select  the  fluid  itself.  For  this  purpose 

(4,5) 

we  have  considered  hydrogen  as  the  primary  working  flviid  with  seedant 

additives  of  three  aUcali  metal  vapors,  Na,  K  and  Cs.  Ihe  particular 

ccnposition  we  have  selected  is  a  mixture  of  85%  H  ,  5%  Na,  5%  K  and  5%  Cs 

2 

by  volume.  Hjdrogen  provides  low  molecular  wei^t  (and  hi^  ^)ecific 

iitpulse)  vhile  the  three  alkali  metals  provide  absorption  capabilities  in 

ccsiplementary  peirts  of  the  soleu:  spectrum.  Absorption  in  alkali  metals  is 
(5) 

conprised  of  electronic  transitions  in  dimers  at  the  lower  temperatures 
(below  about  2500  K)  and  of  photoionization  at  hi^er  temperatures  (2000  to 
3000  K) .  The  characteristic  resonance  doublet  of  alkali  metals  also 
contributes  to  absorption.  Inverse  bremsstrahlung  generally  does  not 
became  significant  until  above  solar  tenperatures. 

Ihe  net  absorption  coefficient  for  the  5%  mixture  of  Na,  K,  and  Cs  in 
tQdrogen  is  given  in  Fig.  1  for  wavelengths  corresponding  to  solar 
radiation.  The  multiple  absorption  paths  in  alkali  meteds  and  the 
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different  wavelength  d^jendencies  of  the  three  alkali  elements  leads  to  a 
very  ocnplicated  absorptivity-tenperature-wavelength  dependency  as  Fig.  1 
shews.  In  peurticular,  the  absorptivity  increases  with  terrperature  at  some 
wavelengths,  but  decreases  with  teirperature  at  others.  Ihis  latter  effect 
arises  becaiase  of  dimer  absorption  and  reflects  the  fact  that  the  number 
density  of  dimers  decreases  with  tenperature  because  of  dissociation.  For 
cenputations,  the  absorptivity  curve  of  Fig.  1  was  divided  into  wavelength 
bands  in  each  of  v4iich  the  absorptivity  was  taken  as  constant.  The  number 
of  bands  used  was  typically  fron  six  to  ten. 

Calculations  of  heat  absorpticxi  in  a  flowing  gas  were  cenputed 

(4,5) 

numerically  using  oontenporary  numerical  methods  .  The  ccitputations 

iise  a  focussed  solar  beam  with  ^X3t  size  chosen  to  correspond  to  a 

concentration  ratio  of  20,000  to  one.  Such  concentration  ratios  correspond 

to  a  reascnable  ufper  limit  on  the  concentration  factors  that  can  be 

achieved  in  practice.  Seme  r^resentative  results  for  flow  in  a  straight 

duct  and  for  flow  in  a  converging  cijct  that  parallels  the  convergence  of 

the  focussed  solar  radiation  are  given  in  Figs.  2-5.  Figures  2  and  3  show 

tenperature  contours  in  the  straight  and  converging  sections,  respectively, 

-4 

for  an  inlet  velocity  corresponding  to  a  Mach  number  of  4  x  10  (about 
0.75  iiv/s) .  For  this  case,  the  dimer  bands  lead  to  rapid  absorption  and 
ste^  tenperature  gradients  near  the  vpstream  windew.  For  the  strai^t 
duct  case,  about  75%  of  the  inccraing  radiation  is  absorbed,  but  a  large 
portion  of  this  is  re-radiated  to  the  walls  so  the  net  fraction  absorbed  is 
only  about  25%.  Much  of  this  re-radiation  to  the  walls  is,  hewever, 
recovered  by  using  the  walls  as  a  regenerative  heat  exchanger.  Ihe 
non-absorbed  radiation  also  goes  to  heat  the  walls  so  it  too  can  be  used 
for  regenerative  purposes.  All  in  cill,  the  radiative  flux  an  the  walls  is 


(1,2,6) 

nucii  less  than  in  the  corresponding  direct  solar  heated  case  and  the 

peak  fluid  tenperature  is  substantially  hi^er.  Seme  fraction  of  the 
re-radiated  energy,  however,  escapes  out  the  solar  window  so  that  the  net 
efficiency  is  estimated  at  about  75%. 

Ihe  converging  geometry  results  in  Fig.  3  show  a  sli^tly  hi^er  peak 
tenperature  of  3700  K  but  the  energy  balances  are  more  or  less  similar  to 
those  for  the  strcii^t  duct  case  in  Fig.  2. 

Ihe  power  cotpled  directly  to  the  plcisma  d^sends  on  the  inceming  flow 

velocity.  To  test  this  sensitivity,  we  cenputed  a  case  for  a  speed  that 

was  twice  that  for  the  results  shown  in  Figs.  2  and  3.  Ihe  results  for 

-4 

this  hi^er  case  (M  =  8  x  10  ,  u  =  1.5  Vs)  are  given  in  Fig.  4.  Ihe 

hi^er  speed  decreases  the  peak  temperature  sli^tly  (to  about  3500  K)  but 
increases  the  net  energy  absorbed  to  about  30%.  This  gives  an  idea  of  the 
type  of  control  that  is  available  by  changing  the  fluid  velocity. 

As  a  final  effect,  we  show  in  Fig.  5  the  effect  of  using  a  solar 
concentration  of  10,000  to  one  rather  than  the  20,000  to  one  ceise  shown 
above.  Ihe  reduced  concentration  ratio  decreases  the  peak  tenperature  from 
3500  to  3000  K.  Ihe  net  energy  absorbed  is  eilso  slii^tly  lower  than  for 
the  hi^er  concentration  ratio  case. 

Ihe  overall  OOTiclusion  of  the  direct  absorption  modeling  is  that  solar 
energy  can  be  absorbed  directly  in  a  flowing  gas  by  using  dimer  seedants  in 
a  primary  hydrogen  flow.  The  resulting  system  gives  much  lower  material 
tenperatures  and  substantially  hi<^er  gas  tenperatures  than  for  the  direct 
absorption  case.  The  dimensions  of  the  absorber  are,  however,  quite  large 
for  the  thrust  sizes  considered.  Solar  roedeets  with  thrust  capabilities  of 
a  few  hundred  Newtons  and  specific  iitpulses  of  500  to  1000  sec  appear 
feasible.  More  detailed  parametric  studies  are  necessary  in  conjunction 
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with  experimental  verification  to  quantify  these  generail  conclusions  and  to 
optimize  the  cycle. 

B.  Analytical  Modeling  of  Microwave  Plasmas 

A  detailed  stucfy  of  microwave  plasmas  has  also  begun  using  a  method 
similar  to  that  developed  for  the  solar  problem.  Ihis  analytical  study  is 
to  augment  the  companion  experiment  of  Micci  described  in  the  first  part  of 
this  r^xjrt.  Ihe  equations  governing  the  interaction  between  a  microwave 
field  cuid  a  flowing  gas  cure  the  compressible  Navier-Stokes  equations  for 
the  fluids  and  the  Maxwell  equations  for  the  electromagnetic  field.  We  are 
presently  considering  the  solution  of  these  equations  in  axisymmetric 
geometries.  The  equations  eure  formulated  in  generalized  bocfy-fitted 
coordinates  that  allow  arbitrary  geometries  and  necessary  grid  stretching. 
Time-marching  methods  are  used  to  solve  the  equations. 

Ihe  configurations  of  interest  center  around  the  floating  plasmas 
described  in  the  first  part  of  the  flew.  To  date,  we  have  detained 
solutions  with  specified  heat  addition  in  the  sphere  cylinder  cenibination 
used  for  the  ei^jeriments  and  coupled  plasma-gas<^T:iamic  solutions  in 
cylindrical  ducts  (see  Figs.  6  and  7) .  To  date,  these  computations  omit 
re-radiation  losses,  but  do  take  into  account  temperature  d^)endent 
viscosity  and  thermal  conductivity.  Effects  of  gravity  have  also  been 
included  (with  a  downward  gas  flow  to  mimic  the  ejperiments) ,  but  the 
results  presented  herein  omit  this  effect.  Figure  7  shows  tenperature 
contours,  streamlines  and  the  specified  heat  addition  contours  for  a 
Reynolds  number  of  400  and  a  heat  flux  of  400  watts. 

Figure  7  shows  the  outline  of  the  micrewave  cavity  throu^  vAiich  the 
fluid  flows  in  a  smaller  quartz  tube.  Ihe  tube  is  here  taken  as 
cylindrical.  Ihe  original  undisturbed  standing  waves  in  the  cavity  for  the 
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axisynroetric  case  are  given  in  Fig.  8.  Hie  corresponding  microwave 
contcurs  in  the  presence  of  a  conAicting  medium  are  given  in  Figs.  7  and  9 
for  two  different  flew  conditions.  Clearly,  the  presence  of  a  conducting 
gas  and  of  energy  addition  completely  alters  the  micrewave  pattern. 
Temperature  contours  in  these  two  figures  shew  that  peak  tenperatures  of 
9600  and  9900  K  etre  reached  for  heat  addition  rates  of  2.8  and  3.7  ]<W.  A 
summary  of  four  cases  cesnputed  to  date  is  given  on  Table  1. 

Current  enphasis  is  on  coiputing  coipled  flcwfields  in  the  exact 
geometry  of  the  esperimental  apparatus  using  the  same  and  epposite 
gravitationcil  orientation.  Specific  problems  being  addressed  include 
methods  for  mecisuring  the  reflected  power  frem  the  cavity  and  methods  for 
modeling  the  line  radiation  loss  from  the  plasma. 
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TABLE  1 


mYSICAL  aMDITICTIS  FOR  KICROWAVE-GASDYNAMIC  CAinJIATIONS 


CASE  1 

CASE  2 

CASE  3 

CASE  4 

Re  # 

90 

450 

180 

18 

m  kg/s 

3  X  lO"^ 

15  X  10“'* 

6.5x10“"^ 

6.5xl0"^ 

Inlet 

velocity 

n/s 

3.9 

3.9 

0.6 

0.06 

Heat 

added,  1<W 

2.8 

3.7 

1.7 

0.27 

Maximum 

Tenperature 

9590 

9890 

9690 

8670 

Total  Pressure 

N/m 

0.2  X  10^ 

1.0  X  10^ 

3.0  X  10^ 

3.0  X  10^ 

T.  =  1000  K;  gas:  helium;  flew  geanetry:  0.1  m  diameter, 

0.2^  m  long;  cavity  gecroetry;  0.18  m  diameter,  0.17  m  long 
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AnSOnPTlON  COEFKICIKNT  ( 


Fig.  1  Complete  absorption  coef 
Cs  in  Hydrogen.  Pq  is  3 


g-  3  Temperature  contours  for  the  converging  duct  case. 
R  =  20,000,  M  =  4  X  10'^. 
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Fig.  5  Temperature  contours  for  the  converging  duct  case 
R  =  10,000,  M  =  4  X  10"^. 
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Fig.  6  Effect  of  specified  heat  addition  on  flow  through  a 
sphere-cylinder  combination.  Re  =  AOO;  Q  =  1500  W. 
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Fig. 


Standing  microwave  pattern  in  cavity  in 
electrically  conducting  medium. 
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emperature  contours  (dashed  lines;  and  microwave  potential  lines  (so 
ines)  in  coupled  mi crowave-gasdynami c  interaction.  Re  =  450,  Q  =  3- 
=  15  X  10“^  kg/s.  (Case  2). 


